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Abstract 

Fluorescence light is induced by extensive air showers while developing in the Earth's atmosphere. The number of 
emitted fluorescence photons depends on the conditions of the air and on the energy deposited by the shower particles 
at every stage of the development. In a previous model calculation, the pressure and temperature dependences of the 
fluorescence yield have been studied on the basis of kinetic gas theory, assuming temperature-independent molecular 
collision cross-sections. In this work we investigate the importance of temperature-dependent collision cross-sections 
and of water vapour quenching on the expected fluorescence yield. The calculations will be applied to simulated air 
showers while using actual atmospheric profiles to estimate the influence on the reconstructed energy of extensive air 
showers. 
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1. Introduction 

Several air shower experiments like HiRes [T], 
the Pierre Auger Observatory [2J, and Telescope 
Array j3], are using the fluorescence technique for 
detecting extensive air showers (EAS) induced by 
ultra-high energy cosmic rays. Measuring the flu- 
orescence light that nitrogen molecules emit after 
being excited by charged particles of EAS is cur- 
rently the most direct method for determining the 
energy of EAS in a model-independent way. A thor- 
ough understanding of the light emission process 
is necessary to obtain the primary energy of EAS 
with high precision. 

In this paper, we extend our previous model cal- 
culation for the fluorescence light emission [4j by 
including the latest results on input parameters and 
their temperature dependence as obtained in labo- 
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ratory measurements. For the reconstruction of air 
shower events, the light emission has to be known in 
dependence on altitude in the Earth's atmosphere 
at which the shower is observed. Up to now, the 
altitude dependence has been considered by includ- 
ing air density profiles and coUisional quenching of 
nitrogen-nitrogen and nitrogen-oxygen molecules as 
described by kinetic gas theory. The cross-sections 
for collisional quenching were taken to be tem- 
perature independent. However, the cross-sections 
describing collisional quenching are known to be 
temperature-dependent [S]. Griin and Schopper [5] 
found a decreasing collisional quenching cross- 
section with increasing temperature. Recently, the 
AirFly experiment has studied collisional quenching 
cross-sections in dependence on temperature [6^. 
These data are also included in the model calcu- 
lations presented in this article. In addition we 
investigate the influence of water vapour on the 
fluorescence yield, using relative humidity measure- 
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ments performed at the site of the Auger detector 
in Argentina. 

2. Model calculation and experimental data 

Following the mathematical description in 1' , the 
fluorescence yield (number of photons of wavelength 
A produced per meter track length) is written as 
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with e\ being the fluorescence efficiency at wave- 
length A without coUisional quenching, p is air pres- 
sure, p'^j, is a reference pressure at which the mean 
life time of the radiative transition to any lower state 
To is equal to that of coUisional quenching Tc- The 
index v' indicates the excitation level of a band sys- 
tem. The air density is given by pair and the energy 
deposited locally by a charged particle of an EAS is 

dE 
dX- 

Assuming air to be a two-component gas, the re- 
lation between p and p'^, is given by 
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with the masses per mole for nitrogen M^.n and 
oxygen Mm,o and the fractional part per volume 
of the two gas components. The temperature depen- 
dence of the coUisional quenching cross-sections is 
parametrised as 
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For example, in case of the AirFly experiment, it 
is given by cr^^ ^, = crj^s^y ■ 293""'^' . The cross- 
sections have been measured for the bands at 
313.6 nm (2P(2-1)), 337.1 nm (2P(0-0)), 353.7 nm 
(2P(l-2)), and 391.4 nm (IN(O-O)). In the calcula- 
tions presented here, it is assumed that ai,' is the 
same for all bands within its band system. In the 
AirFly experiment, a^' have been measured for dry 
air and no differentiation has been made for nitro- 
gen or oxygen. For fitting those data into the model 
calculation, it is assumed that these a-cocfficicnts 
can be applied to the quenching cross-sections of 
both N-N and N-0 collisions. 



An absolute calibration of the AirFly experiment 
has not yet been published. Therefore, the fluores- 
cence efficiency £337. i„m tias been set to 0.082% of 
the deposited energy, the value given by Bunner [7] . 
The same normalisation for £337 inm used in the 
calculations of |4j. 

In the following we will show two model calcula- 
tions based on AirFly measurements. If the temper- 
ature dependence of the coUisional quenching cross- 
sections is not considered, Eq. ^ is used with the 
parameters p'^, given in [8^. Eq. ^ is applied for 
calculating the fluorescence yield with temperature- 
dependent cross-sections. 

The resulting fluorescence yield spectrum, calcu- 
lated for the US Standard Atmosphere at sea level, 
can be seen in Fig. [1] in comparison with calcula- 
tions and other measurements presented in detail 
in [3]. The calculations following the mathemati- 
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Fig. 1. Fluorescence yield spectra of several calculations and 
measurements for 0.85 MeV electrons as exciting particles in 
the US Standard Atmosphere at sea level. The calculations 
following the mathematical description in (4] are labelled 
with BK_name, where BK stands for the initials of the cor- 
responding author of |4] and name indicates the authors of 
the input parameters used. Details can be found in the text. 
The bars indicate the combination of calculation and input 
parameters which is favoured in [4]. 

cal description in [4j are labelled with BK_name, 
where BK stands for the initials of the correspond- 
ing author of and name indicates the authors 
of the input parameters used. Input parameters in 
Eq. [2] are the deactivation constants which are the 
radiative life time tq.u' and the coUisional cross- 
section between nitrogen and nitrogen molecules 
cnn.i/' and between nitrogen and oxygen molecules 
c^NO.u'- Bunner provides coUisional cross-sections 
and radiative life times for the most prominent 
band systems of nitrogen 7^ . Using these input pa- 
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rameters in the calculation, the results are labelled 
with BK_Bunner. Recent measurements by Moro- 
zov et al. [9j were performed for the 2P v' — 0,1 
band systems. For results named BK_Morozov, the 
values from Bunner are replaced by the newer data 
by Morozov et al. where available. An alternative 
calculation of the fluorescence efficiency without 
coUisional quenching e° is also presented in [4] . Here 
the Einstein coefficients A^/^// and the radiative 
life times from Gilmore et al. [lOj and the relative 
apparent excitation cross-section Qapp from Fons 
et al. [TI] are used. For details of this procedure 
see [4]. The resulting fluorescence yield is labelled 
with BK.G.-F. 

For a relative comparison of 19 bands between 300 
and 400 nm see Fig. O The basis is the calculation 
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Fig. 2. Relative comparison of 19 bands of the 
"BK_Morozov"-calculation with measurements and further 
calculations. The absolute fluorescence yield of these contri- 
butions can be seen in Fig. [T] with the same colour labels. 

"BK_Morozov" . 

3. Temperature-dependent cross-section 

In the following, the fluorescence yield for 
0.85 MeV electrons and for EAS are studied us- 
ing the "BK_Morozov" -calculation upgraded by 
temperature-dependent coUisional quenching cross- 
sections. 

The fluorescence yield profiles for a 0.85 MeV 
electron, corresponding to an energy deposit of 
0.1677 GeV per g/cm^, is shown in Fig. [3] for the 
US Standard Atmosphere. The yield is the sum of 
26 bands between 300 and 400 nm as listed in 8 . 
The solid black line represents the results where no 
temperature-dependent coUisional quenching cross- 
sections have been considered. The red dashed line 
includes the temperature-dependent cross-section 



4.5 
4 

3.5 



- 










L AirFly 










- AirFly v 


I'ilh tcmp.-tlL'p. ur 


jss-sctlion 






■ . . . . 1 . . . 


tilh temp. -dep. L:r 
1 .... 1 


... 1 ... . 


n iind S (J hopper 
1 .... 1 . 





height a.s.l. (km) 

Fig. 3. Fluorescence yield profiles for a 0.85 MeV electron 
in the US Standard Atmosphere with measured parameters 
from the AirFly experiment 8 G . The yield is the sum of 26 
bands between 300 and 400 nm as listed in [H]. See text for 
details. 

as measured by AirFly. The blue dotted line reflects 
the AirFly data where the a-coefficients for the 
temperature-dependent cross-sections are replaced 
by the data from Griin and Schopper [5] . The value 
has been extracted from Fig. 6 of [5] and is there- 
fore quite imprecise. Additionally, a caveat has to 
be applied to this comparison: Griin and Schopper 
measured one a-coefficient for the entire wavelength 
range in pure nitrogen and not in air. For this simple 
comparison, the value for N-N collisions has been 
used for the p' of air. However, the AirFly result 
confirms a decreasing cross-section with increasing 
temperature and also the absolute scales of the two 
publications are in quite good agreement. The new 
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Fig. 4. Relative difference of the fluorescence yield calcu- 
lated with temperature-dependent coUisional cross-sections 
compared with the former model calculation. 

dependence reduces the fluorescence yield at sea 
level by about 2% compared with the former model 
calculation, see Fig 21 This deficit increases with 
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increasing altitude to about 10% at 11 km a.s.l. in 
case of AirFly and to about 20% in case of Griin 
and Schopper data. 

In the next step, the calculation of the fluores- 
cence yield is performed using Argentine atmo- 
spheres as given in [12j . The former model calcu- 
lation leads to fluorescence yield profiles for the 
seasonal atmospheres as can be seen in Fig. [5] T . 
Applying the temperature-dependent coUisional 
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Fig. 5. Fluorescence yield profiles for a 0.85 MeV electron 
in the US Standard Atmosphere and measured Argentine 
atmospheres as given in |12) . The given yield is a sum of all 
emitted photons between 300 and 400 nm. 

cross-sections of AirFly in the calculations using 
Argentine atmospheres, the fluorescence yield pro- 
files are strongly distorted, see Fig. [S) The diflFerent 
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Fig. 6. Fluorescence yield profiles for a 0.85 MeV electron in 
the US Standard Atmosphere and measured Argentine at- 
mospheres as given in 1121 . The curves are calculated with 
the "BK_Morozov" -model of [i] combined with the tempera- 
ture-dependent coUisional quenching cross-sections from Air- 
Fly [g. 

temperature profiles lead to a varying strength 
of the influences of the temperature-dependent 
cross-sections. In general, it can be stated that 



the temperature dependence reduces the increase 
of the fluorescence yield within the lowest 11 km 
in the atmosphere significantly. Furthermore, the 
differences between the seasonal atmospheres are 
reduced compared with the calculations without 
the temperature-dependent cross-sections. Fig. [7] 
displays the difference of the fluorescence yield in 
Argentine atmospheres to that in the US Standard 
Atmosphere with the former calculation. 
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Fig. 7. Difference of the fluorescence yield profiles for a 
0.85 MeV electron in measured Argentine atmospheres to 
those in the US Standard Atmosphere without tempera- 
ture-dependent coUisional quenching cross-sections. 

For estimating the importance of temperature- 
dependent cross-sections on reconstructing FAS 
profiles, two average iron- induced EAS with Eq = 
10^9 eV have been simulated with CORSIKA [13], 
one with vertical incidence and the other with 60° 
inclination. The simulations have been performed 
with the US Standard Atmosphere and afterwards, 
the conversion from atmospheric depth X to geo- 
metric altitude h has been done applying Argentine 
atmospheres [12]. For these EAS profiles, the fluo- 
rescence light is calculated and shown in Fig. [8] for 
the vertical case and in Fig. (TU] for the EAS with 
60° inclination. The corresponding differences of 
the fluorescence light in Argentine atmospheres to 
that in the US Standard Atmosphere are displayed 
in Figs. [9] and [TT] The primary energy of EAS can 
be obtained be converting the fluorescence light 
into local energy deposit and then integrating over 
the entire profile of an EAS. Including temperature- 
dependent coUisional quenching cross-section in the 
calculations, the expected shower light profile of an 
EAS with vertical incidence is reduced by 2.7% in 
Argentine summer, by 3.2% in autumn, by 3.7% 
in spring, and by 4.8% in winter. Doing the same 
for the 60° inclined shower, the reduction of the 
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Fig. 8. Fluorescence light profiles for an iron-induced exten- 
sive air shower with Eq = 10^^ eV and vertical incidence in 
the US Standard Atmosphere and measured Argentine at- 
mospheres as given in |12) . The fluorescence emission is cal- 
culated with the "BK_Morozov" -model of [4] combined with 
the temperature-dependent coUisional quenching cross-sec- 
tion from AirFly [6] and the extensive air shower is simulated 
with CORSIKA [13]. 
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Fig. 9. Difference of the fluorescence light profiles as shown 
in Fig. \8\ 

expected light increases to 6.3% in summer, 6.8% 
in autumn, 7.1% in spring, and 7.5% in winter. A 
former study has shown that varying atmospheric 
profiles influence the longitudinal shower develop- 
ment. It yields in slightly distorted profiles of the 
energy deposit of the air shower, which means an 
uncertainty of the energy reconstruction of the pri- 
mary particle of less than 1% [T3]. The position of 
the shower maximum has also been studied and a 
shift of about -15 g cm~^ on average was found. 
Applying additionally the temperature-dependent 
cross-sections, the position of the shower maximum 
is only shifted slightly beyond it. In all atmospheric 
models, the additional shift is less than 50 m. 

The same calculations have been performed for 
proton-induced air showers with the same parame- 



Fig. 10. Fluorescence light profiles for an iron-induced ex- 
tensive air shower with Eg = 10^^ eV and 60° inclination 
angle. Apart from the inclination, everything is the same as 
in Fig. E] 
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Fig. 11. Difference of the fluorescence light profiles as shown 
in Fig. [To] 

ters. The reduction of the expected light is increased 
by about 0.5% compared to the numbers of the iron- 
induced showers. 



4. Vapour quenching 

In this study, the model calculation has been 
expanded by including coUisional quenching due 
to water vapour. An additional term is inserted in 
Eq. ^ to account for the collisions between nitro- 
gen and water vapour molecules. The corresponding 
cross-section has been measured by e.g. Morozov et 
al. [9J and Waldenmaier ,15]. It has to be stressed 
that no temperature-dependence has been mea- 
sured for the coUisional quenching cross-section 
between water vapour and nitrogen. Typically in 
all atmospheric models, the humidity is set to zero. 
For this study, seasonal average profiles of relative 
humidity at the site of the Pierre Auger Observa- 
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tory in Argentina, measured during night times, are 
fitted, see Fig. [HI 
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Fig. 14. Difference of the fluorescence yield profiles in Ar- 
gentine atmospheres as shown in Fig. [m to the US Standard 
Atmosphere. 



Fig. 12. Seasonal average profiles of relative humidity at 
the site of the Pierre Auger Observatory in Argentina as 
measured during night times. 

Firstly, vapour quenching has been considered in 
the former model calculations applying the cross- 
sections measured by Morozov et al. [9] without 
temperature-dependent cross-sections. The result- 
ing fluorescence yield profiles for a 0.85 MeV elec- 
tron can be seen in Fig. 1131 Vapour quenching re- 
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Fig. 13. Fluorescence yield profiles for a 0.85 MeV elec- 
tron in the US Standard Atmosphere and Argentine atmo- 
spheres. The fluorescence emission is calculated with the 
"BK_Morozov"-model of [3] including vapour quenching us- 
ing the humidity profiles given in Fig. 1121 however no tem- 
perature-dependent coUisional quenching cross-section has 
been applied. 

duces the fluorescence yield mainly in the lower part 
of the atmosphere, compare with Fig. [3 because the 
water vapour content in the atmosphere is low at 
higher altitudes. The differences of the calculations 
in Argentine atmospheres to the US Standard Atmo- 
sphere are shown in Fig. [141 The effect is largest dur- 



ing summer, because more water vapour can be con- 
tained in warmer air. The reduction of fluorescence 
yield is most significant near the ground, about 20%, 
and becomes less than 5% above 7 km a.s.l. During 
winter, the effect is of minor importance and only 
visible below about 3 km a.s.l. 

Secondly, vapour quenching has been included in 
the model calculation presented in Sec. [H Again, 
the fluorescence yield is determined for a 0.85 MeV 
electron in Argentine atmospheres, Fig.[T21 The dif- 
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Fig. 15. Fluorescence yield profiles for a 0.85 MeV elec- 
tron in the US Standard Atmosphere and Argentine atmo- 
spheres. The fluorescence emission is calculated with the 
"BK_Morozov" -model of 4 including vapour quenching us- 
ing the humidity profiles given in Fig. 1121 and tempera- 
ture-dependent cross-sections. 

ference of the fluorescence yield in Argentine atmo- 
spheres to that in the US Standard Atmosphere can 
be seen in Fig. 1161 The additional vapour quenching 
changes the profiles shown in Fig. [6] mainly in the 
lowest part of the atmosphere. In summer, the fluo- 
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Fig. 16. Difference of the fluorescence yield profiles as shown 
in Fig. [15] 

rescence yield is reduced significantly and in winter 
the effect is smallest. 

Thirdly, the model calculation including all de- 
pendences is applied to the average iron-induced 
EAS which are used already in Sec.[3l Since the ef- 
fect of water vapour quenching is most important 
near ground, here only the fluorescence light pro- 
files are shown for the vertical shower. Fig. [iTl The 
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Fig. 17. Fluorescence light profiles for the iron-induced ex- 
tensive air shower as shown in Fig. |8] For the fluorescence 
emission calculations, additionally the vapour quenching has 
been included. 

graph of the differences of the fluorescence light in 
Argentine atmospheres to that in the US Standard 
Atmosphere can be seen in Fig. [181 The expected 
light of EAS is reduced by about additional 8.2% 
due to added water vapour quenching in Argentine 
summer, by 5.5% during autumn, by 3.4% during 
spring, and by about 2% during Argentine winter. 
In total, including the temperature-dependent coUi- 
sional cross-sections and the water vapour quench- 
ing, the expected light of the EAS is reduced by 



Fig. 18. Difference of the fluorescence light profiles as shown 
in Fig. [17] 

11.1% during summer, 8.9% during autumn, 7.3% 
during spring, and 6.8% during winter. 

For the 60° inclined shower, the additional effect 
due to vapour quenching is smaller and ranges be- 
tween 1.2% in summer and 0.2% in winter. Com- 
bining the two effects, the expected light is reduced 
by 8.4% during winter, 8.1% during spring and au- 
tumn, and 8.0% during summer. 



5. Conclusion 

The effects of temperature-dependent coUisional 
quenching cross-sections and of quenching due to 
water vapour have been studied. Both effects lead 
to a significant reduction of the fluorescence yield 
in the lower part of the atmosphere. Applying these 
calculations to simulated EAS, a distortion of the 
longitudinal shower development is found. A reduc- 
tion of the emitted ligth is expected, which varies 
from about 7% to 11% depending on seasonal at- 
mospheric model and on zenith angle of the EAS. 
Hence, accounting for these effects in the reconstruc- 
tion of the primary energy of EAS, the primary en- 
ergy will be increased by this amount as compared 
with the former model calculations. The position of 
the shower maximum is hardly shifted, in all atmo- 
spheric models the shift is less than 50 m. 
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